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On an exact criterion for choosing the hopping operator 
in the four-slave-boson approach 
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We consider the JV-component generalization of the four-slave-boson approach to the Hubbard 
model, where 1/iV acts as the small parameter that controls the fluctuations about the saddle point, 
and address the problem of the appropriate choice of the bosonic hopping operator z icr . By suitably 
reorganizing the Fock space, we show that the square-root form for Zi, a (originally introduced by 
Kotliar and Ruckenstein) reproduces the exact independent-fermion ((7 = 0) results not only at the 
mean-field (N = oo) level but also to all orders in the 1/N expansion, provided one relaxes the 
usually adopted normal-ordered prescription for Zi, a . This ensures that Zi, a needs not be modified 
at successive orders in the fluctuation expansion, and implies that all correlation functions are 
correctly recovered in the U = limit, a nontrivial result for the slave-boson approach. In addition, 
it provides a stringent requirement on the form of Zi,&, which might be also generalized to alternative 
slave-boson formalisms (like the spin-rotation- invariant formulation). 
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The four-slave-boson method has been widely used 
in the past few years to deal with the Hubbard model 
(with on-site repulsion U). This method is usually im- 
plemented via a mixed fermion-boson functional inte- 
gral, which allows for a systematic expansion in terms 
of fluctuations about an appropriate mean field. The 
fluctuation expansion is, in turn, controlled by introduc- 
ing an additional fcrmion degeneracy N and using 1/N 
as the expansion parameter .El Already at the mean-field 
level, the method provides a reasonable description of 
the phase diagram and of several static quantities for the 
Hubbard model, as shown by comparisons with Quan- 
tum Monte Carlo calculations.cl Since the slave-boson 
approach should be, by its conceiving, appropriate to 
the strong-coupling regime, it might not be expected to 
reproduce the noninteracting fermion limit at the same 
time. In this respect, Kotliar and Ruckenstein (KR) pro- 
posed choosing the hopping operator Zi^ a (which, due 
to the redundancy of the slave-boson Fock space, is not 
uniquely defined) to obtain the corxect result at the 
mean-field level in the U — > limitS It has, however, 
later been questioned whether the same form of z^ a re- 
mains appropriate when corrections beyond mean field 
are introduced.cl It turns out, in fact, that unphysical re- 
sults occur-when the KR form for is used beyond 



mean field.ffi 

It has been recently shown that the above anomalies 
stem from an inappropriate operator ordering of z^o-.u In 
fact, previous to the work of Ref. [j], all slave-boson calcu- 
lations have been invariably carried out by assuming the 
normal-ordering prescription for Zi, a , since normal order- 
ing allows for a direct functional-integral representation 
of the Hamiltonian. If one adopts, instead, the KR form 
of Zi^c without modifying the original operator ordering, 
fluctuation corrections do not spoil the U = mean-field 



results (at least, for the free energy and related static 
quantities), in contrast to the normal-ordered version for 
Zi M . The question naturally arises whether the results 
of Ref. ^, which have been obtained up to leading order 
beyond mean field, also hold when higher-order correc- 
tions are considered, and not only for static but also for 
dynamic quantities. 

In this paper, by exploiting some exact properties of 
the hopping operator and of the related constraints for 
arbitrary values of the degeneracy parameter N, we show 
that the U = independent-fermion result can be re- 
covered (at zero temperature) to all orders in the 1/N 
expansion, provided the strict square-root form for z ia is 
used instead of its normal-ordered version. This result 
holds not only for the ground-state energy, but also for 
the ground-state expectation values of operators written 
in terms of the physical fermion operators (obtained as 
the product of zi a with a pseudofermion operator). 

We consider the TV-component generalization, of the 
slave-boson single-band Hubbard Hamiltonian,lilu of the 
form 



iV 



cr=±l S=l 

+Uj24di (!) 



supplemented by the constraints 

4 d i+ J2 s L^ + e\e l= N , 



<7=±1 



JV 



fls,ahs,o = s{ a Si^ + d\di 



(2a) 
(2b) 
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which are straightforward generalizations of the N = 1 
case considered by KR. Although for N > 1 the con- 
straints (||) no longer guarantee a one-to-one correspon- 
dence between the original fermion and the fermion- 
boson problems, we consider this not an issue since 1/N 
has been here introduced merely as an expansion param- 
eter and the N = 1 value will eventually be selected. In 
(|l|) and (g) ei, Si^ a , and di are slave-boson operators asso- 
ciated with empty, singly, and doubly occupied states, in 
the order, fi t s,a is a pseudofermion operator with spin a 
and component S, and i {%') extends over all lattice sites. 
For arbitrary N, the presence of the bosonic operator z^ a 
in (Q) is required to leave the subspace identified by the 
constraints (^) invariant; its choice, however, is to some 
extent arbitrary due to the redundancy of the Fock space. 
Kotliar and Ruckenstein .exploited this arbitrariness and 
suggested using the formd 



|n e ,n-|-,nj.,n d > 



(5) 



T Ri. a di 



R 



with 



and 



H,a ■ "i,er 



(3) 



(4a) 



R SQ = 



\J N - d\d l - s\ a Si ta 



N - e\ei - s] Si 



(4b) 



in order to reproduce the exact U = result at the mean- 
field level. The normal ordering : i?f : in ( fia| ) stems by 
requiring a simple mapping of the Hamiltonian (^) into 
the corresponding Action of the functional integral. 

It has, however, been shown in Ref. |l]that improved re- 
sults are obtained beyond mean field by relaxing the nor- 
mal ordering in Eq. ( j4a| ) (whereby the normal-ordering 
requirements of the functional integral are met by suit- 
ably reordering the operators in ( j4b| ) within the 1/N 
expansion). In the following, we shall not rely on the 
functional-integral formulation; rather, we will prove ex- 
act results for arbitrary values of N using the operator 
form j4b|). Specifically, we will show that: (i) the mixed 
fermion-boson Fock space restricted by (|^) can be split 
into subspaces (labeled by a quantum number Ji) that 
remain invariant under the action of Zi^ when the form 
(|l]) is adopted; (ii) a particular subspace can be identi- 
fied (with Ji — for alH), where a one-to-one correspon- 
dence between matrix elements of the original fermion 
and of the fermion-boson Hamiltonians can be estab- 
lished for all N. We will also argue that, when {7 = 0, 
the ground state belongs to the subspace {Ji — 0}. 

We begin by analizing the properties of the operator 
( ^b| ) in the bosonic Fock space at a given site i. We 
specify a generic basis state in this space via the bosonic 
occupation numbers n e ,ni;,ni, and rid associated with 
e^e, sjs-f sjsj,, and d^d, in the order: 



The operator (4b) is diagonal in the representation (|5|) 
(by contrast, its normal-ordered version (Bat) has nontriv- 
ial matrix elements in this representation). This allows 
us to verify the following commutation relations for the 
operator z i f given by (J3j) with R i:t7 = i?f^, when act- 
ing on any state of the form (|^) (we will omit below, for 
simplicity, the site index i and the label SQ): 



(6) 



plus their Hermitian conjugates. 

The "physical" (fermion-boson) subspace is identified 
by the constraints (||). The set of bosonic states asso- 
ciated with a given fermionic configuration is thus de- 
termined via the pair of fermionic occupation numbers 

N i,a- = Es=i fis,<rhs,o- {? = ±1) [a given pair (N h N t ) 
may, on the other hand, be associated with more than 
one fermionic configurations]. The constraints (||) thus 
associate the pair (N^N^) with the bosonic subspace 
spanned by the basis states 

\N h N l ;n d > N 

= \N-N r -N l +n d ,N< f -n d ,N l -n d ,n d > (7) 



-N h N- 



N t ) inte- 



where n d can take 1 + min(N^ , 7Vj_ , -/V - 
ger values ranging within max(0, —N + N^+N^) < n d < 
min(N^ , Ni). The subspaces { | iV-p , N±; n d >n } with dif- 
ferent values of (iVj, N±) are connected by the operators 
z a and z\. Consider, in particular, the application of the 



operators z^ (with z 
\n,0;n d = 0> N : 



and 



r ) on the state 



z+ \n, 0; n d = 0>n = C^\n ± 1, 0;n d = 0> 



JY 



(8) 



where the normalization constants equal unity only 
for the SQ form @.Q Therefore, the operators z^ act 
as creation and destruction operators in the subspace 
spanned by {|n, 0; n d — >n ; n = 0, • ■ ■ , N}. The above 
conclusions hold, as well, if one exchanges up and down 
spins. Exploiting the commutation relations (^|), one can 
then unambiguously define the subspace of the physical 
space spanned by the (N + l) 2 (normalized) states 

|J = 0;JV T) JV,i)^ =(z\) N T(z\) N i\0,0;n d = 0> N (9) 

with < (iV-f, N±) < N . Here J is a new quantum num- 
ber which will be essential for the following arguments. 

Consider, next, the subspace with = iVj_ = 1 
spanned by the two states 1 1 , 1 ; n d = > jv and 1 1, 1 ; n d = 
1 >n (cf- Eq. (^)), since n d takes the values and 1 
in this case. Recall that the state |J = 0; 1, l)jy = 
z|zj|l, l:n d = 0>n belongs to this subspace. The com- 
plement of this state in the subspace we are considering 
can then be found by Schwartz orthogonalization. Let us 
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denote this state by | J = 1; 1, 1)n ■ Similarly to what we 
have done in Eq. (H), from the state | J = 1; 1, 1)jv one 
can construct a whole set of (N — l) 2 companion states 
as follows: 



1 U = i;l,i: 



(10) 



where cj{ T *2 is a normalization constant which is no 
longer unity as it was in Eq. (O). This difference between 
the values of the normalization constants for the cases 
J = and J > 1 is due to the square-root choice ( pb| ) 
and will be important in the following. It can further be 
shown that the states of the two sets {\J = 0; N^ , N^)jy } 
and {| J = 1;N^,N\)n } are mutually orthogonal. 

One can then proceed by induction and identify 
[N/2] + 1 orthogonal subspaces, labeled by the quan- 
tum number J(= 0, 1, • • • , [N/2]) and spanned by the 
(N - 2 J + l) 2 states 



a 



IJ;^,^ = (z\) N i- J (zl) N i- J \J; J, J) N 



(11) 



with J < (iVj , N±) < N — J. Note that successive appli- 
cations of z^ do not change the value of J. 

There remains to evaluate the constants C(j^ ^ ) (= 

N y by symmetry). One can readily show that the 
quantity 



v(J) = 



C, 



(J) 

(n+l,m) 



c, 



(J) 

(n,m) 



— Ar(J;n,m| z^zt \J;n,m)]y (12) 



is eigenvalue of the operator Zjzj associated with the 
eigenstate | J; n, m)jsr and is independent of to. Consider 
then the trace of the operator Zjzj in the subspace with 
given values of N-[ — n and Ni = m, where, for the sake 
of definiteness, we choose m < n < N — m. This sub- 
space is spanned by to -I- 1 states, which we can specify 
alternatively by {|J;rt, m)^ , with < J < to} or by 
{|n, to; rid >n , with < < to}. In the first (diagonal) 
basis the restricted trace takes the form 



tr (z T zj)(„ !m) = ^2 ' 



from which the eigenvalue o4 can be obtained as the 
difference 



a n a) = tr {Z]z\)(n, m ) - tr (Zt4)(n,m-1) 



(14) 



In the second (nondiagonal) basis, on the other hand, the 
restricted trace can be readily evaluated since the opera- 
tor z^zt admits a simple representation. Comparing the 
two results, one obtains: 



(J) _ (n + 1 — J)(N — n — J) 



(N-n)(n + l) 



(15) 



Note that = 1, as anticipated, and that a„' = 
for the upper value n = N — J (except for J = 0, see 
footnote o). In addition, we have 



< atf+V 



< a 



(J) 



< a 



(J=0) 



= 1 



(16) 



for any pair (n,n') and J > 1. 

In conclusion, we have shown that (at any given lat- 
tice site) the constrained bosonic Fock space (for given 
N) can be split into the direct sum of subspaces spanned 
by the states {| J; N h Ni) N , with J = 0, • • • , [N/2]}. 
Within each J subspace, the operators z a and z\ act as 
"lowering" and "raising" operators with normalization 

constants given by ( |l2] ) and (|l^) (C(j j) = 1> by defini- 
tion) and destroy the extremal states with N a = J and 
N a = N—J, respectively. In this respect, the J subspaces 
bear some analogy with the sets of states obtained in el- 
ementary quantum mechanics by coupling two angular 
momenta. 

We now return to the (physical) mixed fermion-boson 
Fock space. The constraints (||) associate every fermionic 
configuration with given values of N-\ and N^, with 
the whole set of bosonic states {| J; N^, -ZVjJjv ',J = 
0, ■ • ■ , min(7V T , N l: N- 7V T , N- N t )} (and not just with a 
single bosonic state, like for N = 1). This implies that a 
generic basis state of the mixed Fock space can be written 
(at any given lattice site) as the product of a fermionic 
and of a bosonic state, as follows: 



\^(N 1 ,N l )> IJ-N^NJn 



(17) 



where 4>(N^, N^) is a generic fermionic configuration con- 
sistent with the pair (A^, iVjJ. Note that the "physical" 
fermion operator z a fs.a does not change the quantum 
number J when applied to (|l7j). 

Concerning further the Fock space for the whole lat- 
tice, this can also be split into distinct subspaces, each 
identified by a set of integers {Ji} with i ranging over all 
lattice sites. Each subspace specified by the set { Ji} is 
invariant under the slave-boson Hamiltonian with [7 = 0, 
since d\di is the only operator of the Hamiltonian (^) that 
modifies the quantum numbers J;. In particular, within 



(13) t ne subspace { J; = 0} for all i the operators z.; iC7 and 



connect states of the type ((17]) with different iVj and N[ , 
without changing their norm. In this subspace, the U = 
slave-boson Hamiltonian is thus equivalent to its purely 
fermionic counterpart also when N > la This finding 
is particularly relevant if the ground state of the Hamil- 
tonian belongs to the subspace {Ji = 0}. In this case, 
the ground-state expectation value of any operator in the 
purely fermionic representation equals the ground-state 
expectation value of the corresponding operator written 
in the slave-boson representation, provided the operator 
itself leaves the subspace {Ji — 0} invariant. This holds, 
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e. g., for the correlation (Green's) functions, when the 
"physical" fermion operator Zi <a fi t s,<T with Ri a given by 
(fit]) is used. Since all physical quantities in the purely 
fcrmionic representation scale simply with integer pow- 
ers of N when U — 0, the above conclusion implies that 
in the mixed fermion-boson representation all corrections 
(over and above the mean-field result) vanish identically 
at any order in 1/N. This result confirms and extends 
on general ground the numerical results obtained in Ref. 
|l| at the leading order in 1/N. 

There remains to prove that the ground state of (|l|) 
with U — belongs to the subspacc {Ji = 0}, that is, 



E o [{J i = 0}}<E [{J' l }} 



with {j;} ± {j, = o} , 

(18) 



where ^[{Jj}] stands for the lowest eigenvalue in the 
subspace specified by the set { Ji}. A property similar to 
( U| ) can be readily proved for a single-particle Hamilto- 
nian with site-dependent hopping, of the form 



N 



H[{Zi,<j}] — Z^j fj.S.gZi^Zi' .o f; 



i' ,S,<j 



(19) 



Ji . J (1 -n ii<T + Jn ii<r ) 



Tk,<r{l - ni l<r ) N N 2 (l-m^)n 



2 

i.a 



(1 - 2n il<r ) Ji(Ni, 



<N LlT >) 



N 2 



+ 0(N- & ) (23) 



where < Ni iCr >= Nfii^ is the self-consistent ground- 
state value of Ni ta for given {Ji}. In this way Eq.jT^) 
is validated for a sufficiently large (albeit finite) value of 
N. In this sense, we have proved order by order in 1/N 
that the slave-boson representation does not modify the 
exact (U = 0) indcpendcnt-fermion resultsJa 

In conclusion, motivated by the encouraging numeri- 
cal results obtained previously by implementing the 1/N 
expansion for the four-slave-boson method correctlyjj in 
this paper we have analyzed in detail the structure of 
the enlarged slave-boson Fock space for N > 1 and iden- 
tified an intrinsic dynamical symmetry associated with a 
novel quantum number (J). We have thus been able to 
tune the four-slave-boson method at U — (where this 
dynamical symmetry holds exactly) by comparison with 
the independent-fcrmion results. Our finding consider- 
ably limits (and possibly eliminates) the arbitrariness on 
the choice of the hopping operator z^ a . 



=±i s=i 



where tn = (by assumption) and Zi^ are real and posi- 
tive numbers. Let Eo[{Zi^}] be the ground-state energy 
associated with a given configuration {Zi tl7 }. Then 



whenever 



E [{ZiA\ <E [{ZU] 



Z it<T > Z< > 



(20) 



(21) 



for all i and cr, which can be shown by realizing that the 
ground-state expectation value 



< dH\{zy\ > _dE [{Z}] 



dZ,. 



dZ t . 



(22) 



is nonpositive for any chosen i and cr. At this point we 
note that the slave-boson Hamiltonian (fil) with {7 = 
can be mapped, within each subspace {Ji}, onto the 
purely fermionic Hamiltonian (|l9| ) by replacing the c- 
number Zi M with the Hermitian operator {a^^ ) 1 ^ 2 [ob- 
tained by entering the, fermionic operator Ni M in the 
place of n in Eq. ([15|)].u Property ( jig ) would thus follow 
from Eqs. ©, ©, and @ if the a { £\ could be re- 
placed by c-numbers. In particular, this is approximately 
correct for large N where one can approximate 
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